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Abstract 

We estimate the quark condensate in one-flavor massless QCD from 
the known value of the gluino condensate in SUSY Yang-Mills the- 
ory using our newly proposed "orientifold" large- N expansion. The 
numerical result for the quark condensate renormalized at the scale 
2 GeV is then given as a function of a s {2 GeV) and of possible cor- 
rections from sub-leading terms. Our value can be compared with the 
quark condensate in (quenched) lattice QCD or with the one extracted 
from the Gell-Mann-Oakes-Renner relation by virtue of non-lattice 
determinations of the quark masses. In both cases we find quite a 
remarkable agreement. 



1 Introduction 



In non-supersymmetric theories, such as QCD, it is very difficult, if pos- 
sible at all, to perform reliable analytic calculations in the non-perturbative 
regime. In the supersymmetric version of the theory the situation is much 
better, due to holomorphy. In particular the exact gluino condensate [1] can 
be evaluated in Af = 1 Super- Yang-Mills. 

In previous publications [2, 3] we suggested a precise way of copying 
non-perturbative results from a supersymmetric theory to a certain non- 
supersymmetric theory named orientifold field theory (due to its realization 
via orientifold type-0 string theory [4]). The orientifold theory is an SU(iV) 
gauge theory coupled to a Dirac fermion in the two-index antisymmetric 
representation, Q + Q Similarly, one can consider [5] a generalized orientifold 
QCD (or QCD OR ), which consists of Nf flavors of Dirac fermions in the 
antisymmetric representation (Nf > 1). 

Our purpose here is to carry out an explicit calculation of the quark 
condensate in one-flavor QCD anticipated in Ref. [3]. Let us briefly recall 
the idea behind such a calculation. Consider three one-parameter families of 
gauge theories, the above-mentioned parameter being N, of their common 
gauge group SU(iV): 

• Pure Yang Mills (YM) theory also known as gluodynamics; 

• QCD F i.e. standard 't Hooft's extension of QCD at arbitrary N (the 
number of quarks in the fundamental plus anti-fundamental represen- 
tation Nf is kept fixed); 

• QCD A , i.e. the SU(iV) gauge theory with Nf Majorana fermions in the 

adjoint representation. 

In Ref. [3] we made a simple observation that, as N increases, the general- 
ized orientifold theory QCD OR interpolates between the three other theories 
above. Indeed, at iV = 2, the fermions of QCD OR are gauge singlets, and, 
therefore, QCD OR reduces to YM. At TV = 3, QCD OR obviously coincides 
with QCD F . Moreover, at N — * oo the bosonic sector of QCD OR goes into 
that of QCD A — a straightforward generalization [5] of the results of Ref. [2]. 
This last observation becomes particularly interesting for one massless fla- 
vor, since in this case the limiting large-iV theory QCD A is nothing but the 
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super symmetric generalization of YM theory — SYM theory also known as 
super symmetric gluodynamics. We will limit our attention to this case here- 
after. 

A consistency check of the above statements can be made by comparing 
the coefficients of the /3 functions of these theories at different values of JV, 
as well as the anomalous dimensions of the corresponding fermion bilinears 
ifjifj. In Table 1 we present this check for the JVj = 1 case under discussion. 
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Table 1. 

We use the standard definition of the coefficients of the (3 function from PDG, 
see Ref. [6], 

da - oaf \ A 2 A 3 , f,\ 

The coefficients can be found in [7], where formulae up to three loops are 
given. The anomalous dimension 7 of the fermion bilinear operators is 
normalized in such a way that 

(»¥)„ = «"*(**)„, (2) 

and n and Q denote the normalization points. For our present purposes we 
can limit ourselves to the two-loop /3-fu nc t;i ons and the one-loop anomalous 
dimensions. We can check easily that the various coefficients of QCD OR go 
smoothly from those of YM (N = 2) through those of QCD F (JV = 3) to 
those of SYM theory at JV -> 00. 

Since many non-perturbative properties of SYM theory are known, the 
large- JV expansion will provide us with information on the non-perturbative 



2 



behavior of one-flavor QCD at N = 3, modulo 1/N ~ 1/3 corrections. The 
planar equivalence method is applicable to a large class of bosonic correlators 
and can be tested, in principle, in lattice calculations. In this letter we will 
concentrate our attention on a one-point function, the quark condensate, 
which has been explicitly computed in SYM theory, and measured on the 
lattice. 



2 Renormalization-group-invariant quantities 
with a smooth large- N limit 

In order to carry out our calculation we need to define, for each theory, 
renormalization-group-invariant (RGI) fermion-bilinear operators and com- 
pare their vacuum expectation values (VEV's) with the appropriate power 
of the corresponding fundamental RGI scale A. Furthermore, we would like 
always to deal with quantities that can be expanded at large N and fixed 
't Hooft (iV-independent) coupling A, 

q 2 N aN 

as a (possibly asymptotic) power series in 1/N. It turns out that this latter 
requirement calls for a slightly unconventional definition of the above quanti- 
ties since at fixed A a non-integer power of g 2 behaves as a non-integer power 
of iV and is, thus, non-analytic at iV = oo. 

The standard (two-loop) definition of the scale parameter A, which follows 
from the conventions of Ref. [6], is 

16vr 2 A /3l//3 ° 2 / 8tt 2 



(For further details see Appendix A.) Because g 2 is multiplied everywhere 
by (3q, Eq. (4) does not suffer from the above-mentioned problem; it defines 
an ^-independent constant of dimension of mass. In what follows, it will be 
convenient to adopt a more general definition, 



A c = ^(c\(riy^ exp (-£ 

V Po 



(5) 
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where the constant c has a finite large- N limit around which it can be ex- 
panded. In the standard definition c = (3q/(2N), cf. Eq. (4). For the time 
being we will keep c as a free parameter, and will discuss the sensitivity of 
our results to the choice of c later. In a similar manner we introduce RGI 
bifermion operators as 1 

(M) £ = iV- 2 (5\(n)y^ <M>, (6) 

where c, like c, has a smooth large-iV limit. Its impact on (^^g will be 
discussed later. Equation (6) will be applied both to the gluino condensate in 
SYM theory and to the quark condensate in the orientifold theory. With the 
above definitions the condensates and A's approach finite limits as N — > oo. 
Moreover, the ratio 

R W = ^ w 

c 

also approaches a finite limit at large N in both theories, and enjoys a smooth 
1/N expansion. 



3 The gluino condensate and the large- N limit 
of the orientifold condensate 

In SYM theory, where a number of exact results were obtained, the gen- 
eral considerations of Sect. 2 simplify considerably. First of all, the expression 
for A in Eq. (5) becomes exact [8] rather than the two-loop approximation, 

Ai ™ = " 3 (db) exp (~m) ■ (8) 

(In this case the standard value of c = 3/2.) Furthermore, the definition (6) 
of the gluino condensate becomes 

(AA>a = ^ (cA(/*)) (\ a ' a K) ■ (9) 

Note that we deal here with the holomorphic part, there is no complex con- 
jugate term in the right-hand side. The reason is that the exact results for 

1 The kinetic term of the fermion fields is canonically normalized, i.e. C = ^> p^!. 
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the gluino condensate are routinely presented in terms of the holomorphic 
condensate, see below. Generally speaking, the VEV (A a ' a A^) is complex. 
We will assume the vacuum angle 9 to vanish. Then one can choose the 
vacuum state in such a way that (A a ' a A^,) is real. We will discuss shortly 
how the gluino condensate defined in this particular way is mapped onto the 
orientifold theory. 

The exact expression for the gluino condensate in SU(iV) supersymmetric 
gluodynamics can be obtained from weak coupling considerations [9]. All 
numerical factors are carefully collected for SU(2) in the review paper [10]. 
A weak coupling calculation for SU(iV) with arbitrary iV was carried out in 
[11]. Note, however, that an unconventional definition of the scale parameter 
A is used in Ref. [11]. One can pass to the conventional definition of A either 
by normalizing the result to the SU(2) case [10] or by analyzing the context 
of Ref. [11]. Both methods give the same result, see Appendix A. When 
expressed in terms of our A SY m it gives, for the ratio defined in Eq. (7) 



This result is exact, there are no 1/N corrections. Note that all existing 
calculations of the gluino condensate were done in the Pauli-Villars (PV) 
regularization scheme. Equation (10), as it is, holds in that scheme. How- 
ever, as we explain in Appendix B, for SUSY gluodynamics with dimensional 
reduction the PV scheme gives the same result as the more currently used 
MS scheme. Thus, all results following from (10) can be viewed as referring 
to the MS scheme. 

Now let us turn to the orientifold theory which is planar equivalent to 
supersymmetric gluodynamics. The first question to ask is the mapping 
of (A a ' a A^) onto There are many ways to establish a proper 

normalization. The simplest way is the comparison of the corresponding 
mass terms. The Dirac fermion \P of the orientifold theory can be replaced 
by two Weyl spinors, and 77^, so that the fermion mass term becomes: 



Rsym(N) 



cc 



(10) 



m^ty = m^r] + h.c. , 



(11) 



while in softly broken SYM the mass term has the form 




(12) 
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Thus, 

-(AA) <-> (£77), or (AA) <->(##) at = 0. (13) 

The same identification is obtained from comparison of the two-point func- 
tions in the scalar and/or pseudoscalar channels in both theories. 
In the orientifold theory 

^= (l ~f| ( ! + ^ =l + O(l/A0, (14) 

P0 1 + 97V 

and 

A) (1 + g^) 

The nonperturbative planar equivalence [2] implies 

Ron(N) = R SYM K(l/N) , K{l/N) = 1 + 0(1/ jV) , (16) 

where the 0(1/N) terms in K reflect deviations from the SYM/OR equiva- 
lence at non-planar level. 

Expressing the result in terms of the conventional fermion bilinear through 
Eqs. (6) and (10) we arrive at 

x K(l/N) c 1 '^/^ c 1 "^) (17) 

where all quantities refer to those of QCD OR , see Eqs (14) and (15). 

This is our final general result. It shows a dependence on the choice of c 
and even more so on c (since its exponent is of order 1/2). Such dependence 
can be absorbed, however, in the definition of the factor K(l/N) modifying 
just the sub-leading terms. We, thus, rewrite (17) in a simpler form 

K(l/N) = l + 0(l/iV). (18) 
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4 Finite- N corrections and numerical results 



The fact that QCD OR goes into YM theory at N = 2 implies the vanishing of 
the fermion condensate at N = 2. In other words we know for sure that the 
function K(l/N) (as well as the previously introduced K{l/N)) must have 
a zero at N = 2. Moreover, arguments can be given that this zero is of the 
first order. Then we can write 

K{\/N) = U-Q K.(l/N), (19) 

where K*(l/N) is supposed to be free from "large" 1/N corrections. Assum- 
ing that if*(l/3) differs from 1 by ±30% at most, and setting N = 3, we 
arrive at the final formula for the quark condensate in one-flavor QCD 



3 / y Wei / 27 \ 

^ = ^(V3) (A(^) ) exp [-^^ ) 



(20) 



where X(/i) = 3a s ((/,)/2ir, see Eq. (3). 

As has been already mentioned, we expect non-planar corrections in K* 
to be in the ballpark ±1/N. If so, three values for If* (1/3), 

tf,(l/3) = {2/3, 1, 4/3} (21) 

give a representative set. The only thing we need now is the value of A(/x). 
Given //, A^g(/i) and Eq. (21) one can get a numerical evaluation of the 
predicted quark condensate in one-flavor QCD. 

The problem is that one-flavor QCD is different both from real QCD, 
with three massless quarks, and from quenched QCD in which lattice mea- 
surements have been recently carried out [12]. In quenched QCD there are 
no quark loops in the running of a s ; thus, it runs steeper than in one-flavor 
QCD. On the other hand, in three-flavor QCD the running of a s is milder 
than in one-flavor QCD. 

To estimate the input value of A^jg(//) we resort to the following pro- 
cedure. First, starting from a s (M T ) = 0.31 (which is close to the world 
average) we determine A^. Then, with this A used as the input, we evolve 
the coupling constant back to 2 GeV according to the one-flavor formula. In 
this way we obtain 

A(2 GeV) =0.115. (22) 
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Then 



(tyty) = -{0.014, 0.021, 0.028} GeV 3 , /i = 2GeV, (23) 

corresponding to three values of in Eq. (21). 

A check exhibiting the sensitivity of our prediction to the value of A(2 GeV) 
is provided by lattice measurements. Using the results of Ref. [13] referring 
to pure Yang-Mills theory one can extract a s (2 GeV) = 0.189. (Here and 
below everything is in MS.) Then, as previously, we find A^, and evolve 
back to 2 GeV according to the one-flavor formula. The result is 

A(2 GeV) = 0.097. (24) 

The estimate (24) is smaller than (22) approximately by one a. This is 
natural since the lattice determinations of a s lie on the low side, within one 
a of the world average. Using Eq. (24) we would get then 

(M) = -{0.05, 0.07, 0.09} GeV 3 , /i = 2GeV, (25) 

Now we have to compare our prediction with an "empiric" value of the 
quark condensate in one-flavor QCD. Chiral perturbation theory allows one 
to determine the quark masses (see e.g. [14]). The Gell-Mann-Oakes-Renner 
(GMOR) relation [15] then implies 

(M) = -0.015 ±0.005 GeV 3 , /i = 2GeV. (26) 

One should remember that the very basis of this derivation, the GMOR re- 
lation, implies three light flavors. One can hope, though, that this particular 
quantity, ( 1 ^ r ^ r ), is not very sensitive to the number of light flavors, although 
it is difficult to assign any uncertainty associated with the Nf dependence. 

Lattice measurements of fifty) were performed [12] in quenched QCD. 
Two methods were used. The first determination was based on the measure- 
ment of the strange quark mass through a fit of the K-meson mass to its 
empiric value. Then (tyty) was extracted from the GMOR relation. The sec- 
ond determination was a direct measurement of (^^J/). Both methods agree 
as far as the central value is concerned, while the uncertainties are much 
larger in the latter method. We quote here the result [12] obtained in the 
first method, 

(tyty) = -0.019 ± 0.004 GeV 3 , /z = 2GeV. (27) 



8 



If, instead, one uses Eq. (47) of Ref. [12] and substitutes there a in physical 
units from Ref. [13], one gets 

(M) = -0.012 ±0.004 GeV 3 , /i = 2GeV. (28) 

In view of the above, it is not unreasonable to assume that the quark con- 
densate in one-flavor QCD lies between these values. Our educated guess 
is 

(MWfi qcd = -0.016 ±0.005 GeV 3 , /i = 2GeV. (29) 

Comparison with Eq. (23) exhibits a significant overlap! Given all un- 
certainties involved in our numerical estimates, both from the side of su- 
persymmetry/planar equivalence and the "empiric" side, we can state with 
satisfaction that the agreement is quite remarkable. 

5 Conclusion 

We started from supersymmetric gluodynamics where powerful methods, 
such as holomorphy, allow one to exactly calculate the gluino condensate. 
We then applied the non-perturbative planar equivalence obtained in [2] in 
conjunction with the orientifold large-iV expansion [3] to predict the value 
of the quark condensate in one-flavor massless QCD, up to subleading 1/N 
corrections. This seemingly first quantitative application of a 1/N expansion 
in D = 4 produces a value for the quark condensate in remarkable agreement 
with the "empirical value." Hopefully, this may pave the way to a whole new 
line of research based on translating a variety of exact results in supersym- 
metric theories to ordinary one-flavor QCD. 

Two extensions of our method look worth being considered: 

• Evaluation of subleading 1/N corrections, for which we were only able 
to give rough estimates here; 

• Extension of our method in the direction of connecting non-supersymmetric 
or M = 1 supersymmetric theories to Af > 2 theories for which even 
more is known. 

Whether or not either one of these developments can be carried out remains 
to be seen. 
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Appendix A. Master formulae 

In this Appendix we present basic formulae which are repeatedly used as an 
input in the bulk of the paper. 

The master formula for the gluino condensate in SU(A r ) supersymmetric 
gluodynamics is 

(X a a \n = -327r 2 M! v lexp(-|^ (A.l) 

= -AN M| v ie" 1/A , (A.2) 
A 

where M PV is the mass of the Pauli-Villars regulator, g 2 is the coupling 
constant at M PV , the gluino field is normalized in such a way that the fermion 
part of the Lagrangian is ( assuming that the vacuum angle 9 = 0) 

i - ■ 

£ferm = A" -DqqA" . (A. 3) 

r 

For numerical comparisons we need to know Aq 11 ^* 1 . This quantity is 
estimated in a number of ways in Sect. 4. We use the standard formala [6] 
for the running gauge coupling constant at two loops, 

The corresponding expression for A is quoted in Eq. (4). 
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Appendix B. The Pauli-Villars vs. MS regular- 
ization schemes 



All existing calculations of the gluino condensate are performed in the Pauli- 
Villars scheme, while all perturabative calculations and experimental deter- 
minations are routinely carried out in the MS scheme. Therefore, for our 
purposes it is necessary to know the relations between the corresponding a's 
or A's. 

The first derivation of this relation can be found in 't Hooft's pioneer- 
ing paper [16], see Sect. 13. Unfortunately, the key expression (13.7) con- 
tained an error which, unfortunately, propagated in part in some reviews, 
e.g. Ref. [17]. It was corrected by Hasenfratz and Hasenfratz [18], see also 
Ref. [19], as well as in a later reprint of Ref. [16] (see [20]). In pure Yang-Mills 
theory 



The difference is entirely due to the fact that the vectorial index /z of the 
gauge connection A® takes D rather than four values in D dimensions. In 

QCD with N f flavors 1/22 must be replaced by (22 - A^f) \ 



In supersymmetric theories the situation is slightly more complicated on 
the one hand, and considerably simpler, on the other. Indeed, dimensional 
regularization per se cannot be used since it breaks the balance between 
the number of the fermionic and bosonic degrees of freedom. For instance, 
in SUSY gluodynamics the standard dimensional regularization would effec- 
tively imply D — 2 bosonic degrees of freedom and 2 fermionic. 

The problem is fixed by using dimensional reduction with the subsequent 
application of the MS procedure. The supersymmetry is maintained because 
even in D ^ 4 dimensions the numbers of the fermionic and bosonic degrees 
of freedom match. 

The most crucial point can be expressed as follows. In the 't Hooft lan- 
guage the difference between the PV and MS schemes comes entirely from the 
non-zero mode parts of the determinants {quantum corrections in the instan- 
ton background). If one uses a more straightforward perturbative language, 
one can split the calculation of the gauge coupling renormalization (by virtue 
of the background field method, in a weak background) in two parts — the 




(B.l) 
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one associated with the magnetic interaction with the background field, and 
the one associated with the charge interaction. It is easy to see that the 
magnetic part produces no difference between PV and MS. The difference is 
entirely due to the charge part, which is in one-to-one correspondence with 
the non-zero mode parts of the determinants in the instanton calculation. 

For non-supersymmetric theories one must carry out a special dedicated 
calculation to analyze the difference between PV and MS. In supersymmetric 
theories the charge-interaction part in the gauge coupling renormalization 
cancels (by the same token all non-zero mode determinants in the instanton 
background cancel). This cancellation is due to the balance between the 
number of fermionic and bosonic degrees of freedom. Thus, dimensional 
reduction plus MS procedure give rise to the same A as the Pauli-Villars 
regularization, 

Apv = Aj^g, SUSY dimensional reduction . (B.2) 
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